Do metals exist in two dimensions? 
A study of many-body localisation in low density electron gases. 
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Using a combination of ground state quantum Monte-Carlo and finite size scaling techniques, 
we perform a systematic study of the effect of Coulomb interaction on the localisation length of a 
disordered two-dimensional electron gas. We find that correlations delocalise the 2D system. In the 
absence of valley degeneracy (as in GaAs heterostructures) , this derealization effect corresponds to 
a finite increase of the localization length. The delocalisation is much more dramatic in the presence 
of valley degeneracy (as in Si MOSFETSs) where the localization length increases drastically. Our 
results suggest that a rather simple mechanism can account for the main features of the metallic be- 
haviour observed in high mobility Si MOSFETs. Our findings support the claim that this behaviour 
is indeed a genuine effect of the presence of electron-electron interactions, yet that the system is not 
a "true" metal in the thermodynamic sense. 
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Since the discovery of Anderson localisation [1] and 
the subsequent scaling theory of localisation [2] in 1979, 
it has been believed that any disorder, no matter how 
weak, would drive a two dimensional (2D) system to- 
ward an insulating state. This absence of metals in 2D 
has been challenged by the observation [3(] in 1994 of a 
Metal-Insulator Transition in high mobility Silicon MOS- 
FETs. This first set of experiments has been followed by 
similar observations [4j, [5| in a wide range of 2D systems. 
However, as there seemed to be no room for this new 
metal within the widely accepted theory [6j, the experi- 
mental data remained a puzzle to the community. Many 
theoretical scenarios were proposed, rangin g fro m rather 
extreme [TJ, la, |9j to fairly conservative [lfj, [llj . From a 
theoretical point of view, the difficulty in dealing with 
these low density 2D systems is double, as both disor- 
der and interaction must be treated in a non pertubative 
way. Indeed, despite important efforts, even the non- 
interacting problem (Anderson localization in 2D) has 
resisted all analytical approaches so far. 

Here, we take a different route, similar in spirit to what 
has been done numerically [12| to verify the assumptions 
of (non-interacting) scaling theory of localization. Build- 
ing on a numerical technique [l3| that we have devel- 
opped recently, we calculate the many-body localisation 
length of the interacting problem. 

We consider the many-body generalisation of the An- 
derson model used to study localisation. The system 
(with N particles in a lattice made of L x x L y sites) 
is parametrized by two dimensionless numbers that are 
denoted traditionally r s (strength of the Coulomb inter- 
actions) and 1/kpl (strength of the disorder). When the 
electronic density n s decreases, both interaction r s — 
m*e 2 J ' (h 2 e^/im s ) and disorder 1/kpl increase in a fixed 
ratio r\ — r s \fkpl which characterizes a given sample (m* 
effective mass, e electron charge, e dielectric constant, kp 
Fermi momentum and I mean free path) . We distinguish 



three types of systems. Polarised electrons are effectively 
spinlcss and the corresponding system is referred as the 
one component plasma (1-C, M — 1). Non polarised 
electrons (like electrons in GaAs heterostructures) cor- 
respond to the two components plasma (2-C, M = 2). 
The non polarised system with valley degeneracy (elec- 
trons in Si MOSFETs) is the four components plasma 
(4-C, M — 4). The Hamiltonian for the M components 
plasma reads, 
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where el et cp a are the usual creation and annihila- 
tion operators of one electron on site r with inner de- 
gree of freedom er, the sum J2iff')a ^ s restricted to near- 
est neighbours and n? a = cl c? a is the density opera- 
tor. The internal degree of freedom a = 1 . . . M cor- 
responds to the spin (M = 2 for semi-conductor het- 
erostructures like GaAs/AlGaAs) or the spin and valley 
degeneracy (M = 4 for Si MOSFETs). The disorder po- 
tential v? is uniformly distributed inside [-W/2, W/2]. t 
is the hopping parameter and U is the effective strength 
of the two body interaction Vp. To reduce finite size 
effects, Vp, whose expression can be found in 
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obtained from the bare Coulomb interaction using the 
Ewald summation technique. We work at small filling 
factor v = N/(L x L y ) <C 1, where we recover the con- 
tinuum limit. The two dimensionless parameters read 
r s = U/(2ty/mJ) and k F l = 192vw£ 2 /(MW 2 ). 

Without interaction (r s = 0), scaling theory of lo- 
calisation 2] predicts that the scaling function /3(g) — 
d\ogg/d\ogL (which indicates how the conductance g of 
the system changes as one increases its size L) is a func- 
tion of g only. In two dimensions, (3(g) is always negative 
(for arbitrary disorder 1/kpl > 0) so that g always ex- 
trapolates to zero in the thermodynamic limit. It was 




FIG. 1: (3(g) as a function of log g for the 2-C (JV=32, left) and 
4-C system (iV=64, right). The red line is the non-interacting 
scaling function. Symbols stand for r s —0 (triangles), r s =2 
(circles) and r s —4 (squares), for 0.55 < kpl < 24.5. The 
black arrow in the left panel follows a given value of disorder 
for increasing r 3 . 
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FIG. 2: a/i as a function of N for 1-C (left), 2-C (middle) 
and 4-C systems (right panel). The symbols correspond to 
different strengths of interaction (r s = (circles), r s = 0.25 
(down triangles), r 3 — 2 (diamonds), r s = 4 (up triangles) 
and r 3 = 6 (squares)) and of disorder (1 < kpl < 6 for 1-C, 
Kk F l< 2.4 for 2-C and k F l = 1.05 for 4-C). 



found in Ref. [13(] (to which we refer for technical details 
about the numerical method) that, for the 1-C plasma, 
not only the non-interacting [3(g) function is unaffected 
by the interaction, but upon increasing r s , the system 
becomes even more insulating. In this letter, we show 
that the 2-C and 4-C plasma have an opposite behaviour 
from 1-C and get delocalised by electronic correlations. 

Let us start with the study of the (3(g) function, shown 
in Fig. [TJ We do not assume that (3 is a function of g only, 
but we merely plot simultaneously (3 — dlogg/dlogL 
and g while varying both kpl and r s . Upon switching 
on the interaction in the 2-C system, g and (3 first in- 
crease quickly. They reach a maximum at r s rs 2 and 
then decrease. We observe no visible deviation from 
one-parameter scaling: even though g and f3 increase 
with r s , the trajectory [g(r s ),/3(r s )] remains on the non- 
interacting 0(g) curve. In particular, at r s rs 2 where the 
delocalisation effect is maximum, the value of (3 always 
remains strictly negative (3 < so that, even though the 
localisation length £ has increased significantly, it is still 
finite, and the system would still be insulating in the 
thermodynamic limit. The situation for the 4-C system 
is qualitatively different from 2-C as the delocalisation ef- 
fect is much more pronounced: (3 quickly increases with 
r s and reaches = i.e. a metallic state up to our sta- 
tistical precision. While the non interacting points are 
spread over a wide range of g and (3, once the interaction 
is switched on, all the points (but one obtained for the 
strongest disorder) move toward (3 — 0. Hence, the lo- 
calisation length £ becomes much larger than the system 
size and the system becomes practically a metal. Yet, 
as in 2-C, we observe no visible deviation from one pa- 
rameter scaling so that even though £ has increased dra- 
matically we cannot conclude that the system becomes 
a "true" metal in the thermodynamic sense (i.e. a di- 
vergence of £). We note that indications that spin and 
valley degeneracy could play a positive role for transport 
properties were found in the limit of low disorder and 
interaction (kpl — ► oo, r s <C 1) as early as 1983 [laJlq. 



However, these works predicted deviations from one pa- 
rameter scaling which we did not observe in our (non- 
perturbative) regime of parameters (0.5 < kpl < 4 , 
<r s < 10). 

As we observed no deviation to one parameter scal- 
ing, it means that all the relevant information of the 
system is encapsuled in the (one parameter) localiza- 
tion length £(r s ,kpl) of the system. Integrating the 
(3(g) = dlogg/dlogL function, one parameter scaling 
implies that 
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is independent of the size of the system. We have 
carefully tabulated (3(g) without interaction and we use 
Eq.© to extract £. Note that even if there are devia- 
tions from the non-interacting (3(g) , those are small in 
the regime of parameters that we have considered, and 
can be safely ignored for our purpose. We also emphasize 
that this method allow us to extract £ even when £ rs L 
and that we recover g = goe~ L ^ for L>(. In Fig. [21 we 
plot a/t; as a function of N for various values of disorder 
and interaction (a = 1/ 'y/nv is the averaged distance be- 
tween particles). Except for very small N (not shown), 
we find that a/t; is indeed independent of N for all three 
systems 1-C, 2-C and 4-C. Note that for 4-C, we show 
only one value of kpl for different r s = 0, r s = 0.25 and 
r s = 2. For the latter, a/£ = up to our statistical res- 
olution while the corresponding non-interacting system 
is localized. Hence, the observed (strong) delocalization 
effect is robust against increasing N. 

The corresponding localization lengths have been col- 
lected in Fig. [3] which is the chief result of this letter. It 
is, to the best of our knowledge, the first calculation of 
the localisation length in presence of many-body correla- 
tions. 

In the rest of this letter, we turn to the implications of 
these results to the transport properties of 2D gases. A 
strong debate has taken place in the literature on the na- 
ture of the observed metallic behaviour [3j (i.e. decrease 




FIG. 3: Inverse localisation length a/£(r s ) for the 2-C (left) 
and 4-C systems (right) at various N and UfI- For 2-C: 
top to bottom: fc F J=1.05 (iV=32 left triangles, 7V=24 right 
triangles), k F l=1.51 (iV=32 diamonds, iV=24 circles, iV=18 
squares), fci^=2.36 (iV=32 up triangles, Af=18 down trian- 
gles). For 4-C: k F l=0.77 (iV=36 circles), k F l = 1.05 (iV=84 
cross, N=64 down triangles, N=36 up triangles), fcirZ=1.51 
(iV=64 squares, iV=36 diamonds). Inset: zoom of the main 
figure. Lines are guides to the eye. Error bars are roughly 
equal to symbols sizes. 



of resistivity as temperature is lowered). We will argue 
that the above findings provide a simple mechanism that 
accounts for the main experimental observations. The 
first observation of importance in those low density sys- 
tems, is that the Fermi energy is extremely low, so that 
the actual typical temperature is not much smaller than 
the Fermi temperature. This is to be contrasted with 
more conventional systems where only a small fraction of 
the electrons around the Fermi surface take part in the 
transport properties. Hence, the nature of the (highly) 
excited states and in particular their localisation proper- 
ties is of crucial importance here. In a non-interacting 
system, the localization length varies very quickly with 
kpl as electrons gather kinetic energy, typically exponen- 
tially [6j. An immediate consequence is that the excited 
states, with higher kinetic energy, are less localized than 
the ground state. In particular, the polarized system is 
less localized than the (non polarized) ground state. In 
such a situation, one naturally expects a negative in plane 
magnetoresistance: an in plane magnetic field H (with 
no orbital but only Zeeman coupling to the electrons) 
polarizes the system hence delocalizes the system whose 
resistance decreases. Similarly, upon increasing the tem- 
perature T, the less localized excited states get populated 
and the overall resistance decreases. 

Now, we find that this natural order can be reversed 
in presence of interactions. One ends up in a situation 
where the ground state is delocalized while its excited 
states are still strongly localized. In this situation, the 
behaviour of the resistivity upon increasing an in plane 
magnetic field H or the temperature is reversed with re- 
spect to the non-interacting case: populating the highly 
localized excited states leads to an increase of resistiv- 
ity, hence a positive in plane magnetoresistance as well 
as an increase of resistance upon increasing temperature 
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FIG. 4: 5(r 3 , 1/VW) for the 2-C (left) and 4-C system (right 
panel). 8 > (red, +) indicates that the polarized system is 
less localized than the non-polarized one, 8 < (blue, — ) the 
opposite. No data are available in the black regions. 



(as in the "metallic" phase of the Si MOSFETs). Both 
effects are expected to take place when the corresponding 
Zeeman energy g/ig.ff/2 and temperature kT are of the 
order of the characteristic energy on which £ varies, i.e. 
the energy Ep(kpl,r s ) to polarize the system. 

In the "phase diagram" Fig. |4j we have plotted the 
difference S = a/S,NP — a/£,p of the inverse localization 
length for the non-polarized £tfp and polarized £p sys- 
tems. The polarized system is an excited state of the 
non polarized one, and we have verified numerically that 
the picture that follows also apply to partially polarized 
systems (excited states of lower energy). The (r s ,kpl) 
plane can be divided into three regions. In the white 
region kpl > 4, both localization lengths £p and £/vp 
are extremely large so that 8 is essentially zero. Upon 
increasing disorder kpl < 4 at small r s , the localization 
length of the system decreases quickly. Yet, as the inter- 
action is weak, the polarized system is less localized than 
the non-polarized one which leads to a positive S (red, 
+). In the last region (blue, — ), the interaction have de- 
localized the non-polarized system much more than the 
polarized one, hence a negative 6. In our view, this (new) 
blue region corresponds to the regime where the metal- 
lic behaviour has been observed. Note that there is an 
important difference between 2-C and 4-C in the blue 
region: for 4-C, the ground state is (for practical pur- 
poses) delocalized while for 2-C it is not. Hence, for 2-C, 
temperature has two conflicting effects: on one hand it 
activates transport mechanisms such as variable range 
hopping [17J, while on the other hand it populates the 
(more localized) excited states. 

Upon varying density, a given sample moves on Fig. [4] 
according to r s = r\j\Jkpl (straight lines). In the ab- 
sence of interaction, r\ = 2 y /Jle' i / 2 m* /(H 3 ^ 2 eyM), so that 
the mobility /i obtained at low temperature and large 
density (small interaction) can be used as an estimate 
of r\. We note however that upon varying density, r\ 
remains constant only for white disorder. Real disor- 
der is correlated so that r\ is expected to increase at 
low density as larger scales are probed. Also, the dis- 
order seen by the electrons at the Fermi level is partially 



screened by the rest so that the above formula overesti- 
mates r\. Rough estimates are rj rs 10 (high mobility Si 
MOSFET, fx > 10 4 cm 2 .V- 1 .s" 1 ), r) w 50 (electrons in 
GaAs/GaAlAs heterostructures, p, > 10 6 cm 2 .V _1 .s _1 ), 
and rj ~ 200 (holes in GaAs/GaAlAs heterostructures, 
p rs 10 6 cm 2 .V -1 .s -1 ) so that high mobility Si MOSFETs 
are dirtier than their III-V counterparts. 

For definiteness, let us concentrate on the data pre- 
sented in Fig.l of Ref. [10| (hereafter referred as Fig A) 
for a typical high mobility Si MOSFET. At very high 
density (kpl 3> 1, r s <C 1), the system can be described 
within diagrammatic theory which predicts small quan- 
tum corrections to transport properties [6(. This is the 
left part of the white region in Fig. |4] where one ob- 
serves, for instance, the usual signatures of weak local- 
ization. The resistivity p(T) as a function of tempera- 
ture is roughly flat in this region (except at large tem- 
perature where phonons set in). Upon decreasing the 
density, one enters the blue region of Fig. [H at kpl rs 4 
(n s rs 2.5 10 n cm -2 for the sample of Fig. A). In this 
region, one observes the strong positive variation of p(T) 
which has been puzzling the community for more than 
a decade. We attribute this behaviour to the fact that 
in this region, the ground state is (for all practical pur- 
poses) delocalized (see right panel of Fig. [3]) while its 
excitations are not. Hence upon increasing temperature 
those localized states get populated and the resistivity in- 
creases strongly. This scenario implies in particular that 
(i) The corresponding characteristic energy scale where 
the increase of p(T) is expected is the energy Ep(kpl, r s ) 
to polarize the system. We have calculated numeri- 
cally this energy scale in the relevant regime and found 
Ep rs 0.15-Ef where Ep is the Fermi energy. Indeed, one 
observes in Fig. A that T/Tp rs 0.15 corresponds to the 
crossover temperature where the strong increase of p(T) 
takes place, (ii) An in plane magnetic held polarizes the 
system which becomes localized, hence a positive mag- 
netoresistance is expected with a variation of resistance 
of the same order of magnitude as the variation due to 
the change in temperature. This is indeed observed for 
magnetic fields H that correspond to a Zeeman energy of 
the order of Ep (i.e. rs 1.6T at n s = 10 n cm -2 ) as can 
be seen, for instance, in Fig. 10 of Ref. [J]. 

To conclude, we have reported on the first quantitative 
calculation of the localization length in an interacting two 
dimensional electronic gas. The scenario that naturally 
emerges from the data captures all the essential features 
of the experiments at a semi-quantitative level, with no 
adjustable parameters. In particular, our data account 
for: (i) why the metallic behaviour is found at low den- 
sity, i.e. when the corresponding polarized system lies in 
the vicinity of the quantum of conductance, (ii) Why it 
is destroyed by an in plane magnetic field, (iii) Why Si 
MOSFETs, with valley degeneracy, show a much stronger 
metallic behaviour than other materials, (iv) Why disor- 
der matters, and samples with r\ < 1 will not exhibit the 



metallic behaviour (red part of Fig. [4]) while too clean 
samples (white part) will only show a weak signal, (v) 
The characteristic magnetic field and temperature for the 
dependance of the resistivity. 

We have not discussed the metal-insulator transition 
per se, but rather focused on the existence of the metallic 
behaviour. Indeed, for large interaction r s and disorder 
1/kpl, there is no question that the system will even- 
tually become insulating. Although the nature of this 
insulator and of the corresponding transition is clearly 
an interesting topic, the real mystery lies in the nature of 
the metal itself. Our results suggest that this metallic be- 
haviour is due to a non-perturbative interaction delocal- 
ization effect. On the other hand, we found no significant 
deviation from the non-interacting one-parameter scal- 
ing which suggests that the 4-C systems are not "true" 
thermodynamic metals and would become insulating at 
ultrasmall temperature (2-C systems are unambiguously 
insulators). In any case, we find that the question of the 
existence of a "true" metal is not relevant in those ex- 
periments as a different physics (of much higher energy) 
has been probed in practice. 
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